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fN| ' To confront relativity theory with observation, it is necessary to split spacetime 

>' 

into its temporal and spatial components. The (1+3) timelike threading approach 

in ; 

involves restrictions on the gravitational potentials (guv)) while the (3+1) spacelike 
slicing approach involves restrictions on (g^). These latter coordinate conditions 

o 

protect chronology within any such coordinate patch. While the threading coordi- 
nate conditions can be naturally integrated into the structure of Lorentzian geometry 



^ , and constitute the standard coordinate conditions in general relativity, this circum- 

?H ■ 

stance does not extend to the slicing coordinate conditions. We explore the influence 

of chronology violation on wave motion. In particular, we consider the propagation 
of radiation parallel to the rotation axis of stationary Godel-type universes charac- 
terized by parameters r/ > and A > such that for rj < 1 (rj > 1) chronology 
is protected (violated). We show that in the WKB approximation such waves can 
freely propagate only when chronology is protected. 
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I. INTRODUCTION 

Though space and time refer to essentially different aspects of our experience, their unifi- 
cation in the theory of relativity has been a distinct achievement. To interpret observations 
in accordance with this theory, however, we need to split spacetime into its components [l[ . 
Observers are macrophysical entities; therefore, their associated temporal coordinate is ther- 
modynamic time derived from the concept of entropy. This leads to the arrow of time, its 
one-way property that has no analog in space. In principle, classical gravity can change 
this one-way character of time and produce a closed timelike curve ("CTC"). That is, such 
a possibility is not ruled out by the geometric structure of classical general relativity. In 
this paper, we approach the gravitational violation of chronology from the viewpoint of the 
theory of measurement of space and time. 

We begin our discussion of spacetime splitting with ideal inertial observers in Minkowski 
spacetime. Imagine a global inertial frame with coordinates x M = (ct, x) and spacetime in- 
terval —ds 2 = rj^dx^dx", where [f]^) = diag (—1, 1, 1, 1). We define fundamental observers 
to be those at rest in space. The world lines of the fundamental observers are therefore 
temporal coordinate lines that thread spacetime. At each instant of coordinate time t, the 
corresponding orthogonal subspace is the "space" of the static observer. This is schemati- 
cally depicted in panel (a) of Figure 1. This "1+3" splitting is the threading approach to 
separating time and space of a fundamental observer. 

Alternatively, consider any congruence of inertial observers. In their bundle of world lines, 
we identify "local" t = constant hypersurfaces, which must represent the space experienced 
by these observers. The vector normal to such a hypersurface must indicate the direction of 
increasing time coordinate. This "3+1" splitting is the slicing approach to separating space 
and time of an observer family and is schematically depicted in panel (b) of Figure 1. 

To preserve the temporal order of events under Lorentz transformations, we must ex- 
plicitly exclude superluminal phenomena. To see this, consider the temporal Lorentz trans- 
formation from the global inertial frame to another inertial frame moving with velocity v, 
i.e., 

t' = 7(t-^v-x). (1) 

Suppose two events (ti,xi) and (£21 x 2) are causally related; hence, |Ax| < c |At|, where 
Ax = x 2 - x x and At = t 2 - h . Thus, At' = 7 (At - v ■ Ax/c 2 ) implies At'/ At > 0, 



and if t% > t±, then t' 2 > t[. In this way, causality is preserved if there is no superluminal 
propagation. 

The treatment of spacetime splitting can be extended to accelerated observers in 
Minkowski spacetime using the hypothesis of locality, namely, the assumption that an accel- 
erated observer is instantaneously equivalent to an otherwise identical momentarily comov- 
ing inertial observer. Further extension to gravitational fields is then possible via Einstein's 
principle of equivalence. The threading and slicing approaches coincide for inertial observers 
in Minkowski spacetime. Using these extensions to accelerated systems and gravitational 
fields, we will show that they are compatible with each other when a universal temporal 
coordinate exists, but are in general incompatible in a gravitational field that has no cosmic 
time. 

In a gravitational field, the threading and slicing approaches lead to restrictions on coordi- 
nate systems that may be used to cover the spacetime manifold. The threading coordinate 
conditions can be invariantly included in the geometric structure of Einstein's theory, so 
that all spacetime coordinate systems must obey these standard threading admissibility 
conditions. The situation is different with the slicing coordinate conditions; in general, 
they simply impose additional restrictions on coordinate systems so as to exclude causality 
violation via closed timelike curves. 

It is interesting to study the influence of chronology violation on wave propagation. In 
fact, Huygens' principle assumes at the outset that coordinate time monotonically increases 
along the wavefront. Therefore, wave motion may exhibit some unusual features in chronol- 
ogy violating universes. A complete treatment of this topic is beyond the scope of this 
paper; however, we explore here some of the consequences of the absence of a cosmic time 
coordinate for wave motion in a gravitational field. In particular, we study wave propagation 
parallel to the axis of rotation of stationary Godel-type universes. 

Geometrical units (c = G = 1) will be used hereafter. The plan of this paper is as 
follows: In Section II, spacetime splitting is discussed in a gravitational field. The resulting 
conditions on the admissibility of coordinates are examined in Section III. Some aspects of 
wave motion in certain Godel-type universe models are treated in Sections IV - VI. Section 
VII contains a discussion of our results. 



II. OBSERVERS IN A GRAVITATIONAL FIELD: THREADING AND SLICING 

APPROACHES 

A. Threading Approach 

The fundamental observers in a gravitational field are those at rest in space. These static 
observers are naturally connected to a coordinate system. All other observers are pointwise 
related to the fundamental observers by Lorentz transformations. Imagine a fundamental 
observer with four-velocity u in an arbitrary gravitational field, as depicted in Figure EJ 
In general, such an observer is accelerated. The spacetime interval along its trajectory — 
namely, where dx % = 0, for i = 1,2, 3 — denotes its proper time, — ds 2 = gttdt 2 . Therefore, 
we must have g u < 0. Moreover, its orthogonal subspace a represents its local rest space 
and is such that u^dx^ = 0, where u^ = {—gtt)~ 1 ^ >J 'o- Thus the equation that determines 
a is g tt dt + gudx 1 = 0. In general, one can write 

g^dx" ® dx v = g tt (dt + ^dx l ) <g> fdt + ^dx j ) + (g i3 - ^-\ dx i <g> dx j , (2) 

so that 

g^vdx^ (g) dx u | = "fijdx 1 <g> dx\ (3) 

where 

9ti9tj i . \ 

Hi = 9n ~ —+ ■ (4) 

gtt 

It follows from the inverse relationship between (g^u) and (g^) that (7^) and (g^) are 
inverse of each other; similarly, (g^) and (7 U ) are also inverse of each other, where 

f j = 9 ij - ^ • (5) 

The threading approach thus leads to the following coordinate admissibility conditions 

g tt < , (jij) = positive definite matrix. (6) 

Given a system of coordinates that uniquely identify events in a spacetime region, the 
standard threading admissibility conditions have to do with the possibility of existence of 
hypothetical static observers at all events in that region. 

We recall that a real symmetric matrix A is positive definite if and only if there exists a 
real invertible matrix V such that A = V T V. Thus if A is a positive definite matrix, then so 



is its inverse; of course, this can also be seen immediately from the positivity of all of their 
eigenvalues. It follows, for instance, that for a positive definite matrix (Aij), $ = Aijdx l dx^ 
can be expressed in matrix notation as $ = (Vdx) T (Vdx), which is manifestly positive, i.e., 
$ > 0. 

The measurement of spatial distance between a fundamental observer and an arbitrary 
nearby observer by means o 



light signals, as well as the synchronization of their clocks, 
indeed, it is demonstrated there that the element of spatial 
distance, d£, as calculated here, d£ 2 = ^ijdx % dx\ coincides with the radar distance 3j. 



has been discussed in Ref. 



3; 



B. Slicing Approach 

Imagine an arbitrary congruence of observers in a gravitational field, as depicted in Figure 
|3j This bundle of future directed timelike curves could, for instance, represent the point 
particles that form a compact object. The "local" t = constant hypersurface £ is defined 
by dt = 0; therefore, 

g^vdx^ ® dx u \^ = gijdx 1 ® dx J . (7) 

If E is to represent space, (#y) must be a positive definite matrix. Accordingly, its inverse, 
( A f lJ ), must also be a positive definite matrix. Moreover, the normal to S, A/^ = 5^°, must be 
timelike; hence, g^M^Mp < 0, or g tt < 0. Thus the slicing approach leads to the following 
coordinate admissibility conditions 

g tt < 0, (j 13 ) = positive definite matrix. (8) 

These conditions ensure that in the coordinate patch under consideration coordinate time 
increases monotonically along any timelike curve. In particular, the application of these 
conditions to a bundle of fundamental observers implies that in such coordinate systems, t 
monotonically increases along the world line of a fundamental observer. Coordinate systems 
that are admissible according to conditions (IE]) therefore exclude closed timelike curves. 

Let us note that for g u ^ 0, we have in general 

g^d, ® d v = g* (a + £ft) ® (dt + ^a,-) + (> - ££) d z ® dj , (9) 

which should be compared with Eq. (J2J); in fact, a formal duality between them can be 
illustrated by means of the lapse and shift functions. 



space 



t =constant 
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FIG. 1: (a) Schematic plot representing the threading approach to spacetime splitting in a global 
inertial frame in Minkowski spacetime. (b) Schematic plot representing the slicing approach for 
the congruence x = xq + vt, where xq varies over the congruence and v is a constant speed. 




FIG. 2: Schematic diagram representing the future directed world line of a fundamental observer. 
To every admissible coordinate system, one can associate a set of fundamental observers, i.e., those 
at rest in space. The world line of such an observer is a time line. The infinitesimal orthogonal 
hypersurface a is not in general a simultaneity hypersurface for the u congruence. The time lines 
thread the spacetime; hence, this is the threading approach to splitting spacetime into time plus 
space ("1+3"). 




onstant 



FIG. 3: Consider any congruence of future directed timelike world lines V in a standard admissible 
coordinate system. The infinitesimal constant-time hyper surf aces S should all be spacelike. Thus 
their normal vector field should be timelike. The hypersurfaces slice up this spacetime region; 
therefore, this is the slicing approach to the splitting of spacetime into space plus time ("3+1"). 
Note that n^V^ > and the temporal coordinate monotonically increases along such a congruence. 

C. Duality 



In the threading approach, the lapse and shift functions are respectively M and Mj given 



by 



fjti 



M = ^g^ t , Mi = -»—. (10) 

9tt 



Then, we have for metric (J2J) the equivalent representation 

- ds 2 = -M 2 (dt - Midx 1 ) (dt - Mjdx j ) + j ij dx i dx j , 7^ = g tj + M 2 MiMj . (11) 



The form of metric ( TTT1) is adapted to static observers with 4-velocity 

u = ^ , u b = -M(dt - Midx 1 ) , 



(12) 



where the "flat" symbol (b) denotes the fully covariant form of a tensor. In fact, 

— ds 2 = —(u a dx a ) 2 + r ) i jdx % dx^ . (13) 

In a similar way, one may introduce the lapse - shift notation in the slicing approach, 
namely, 



N 



1 

7=F 



N l = - 



9_ 
9* 



(14) 



Normalizing A/" M , so that it becomes a vector of unit length n M , we have 

n = ~{dt - N%) , n b = Ndt . (15) 

Thus in the slicing approach, Eq. ([9]) can be expressed in terms of n a and 7' J in a manner 
that is completely analogous to the threading approach, namely, 

-(d s ) 2 = _(n a a ) 2 + 7*(0i)(0i). (16) 

However, a closer examination reveals that this formal duality does not extend to the physical 
conditions contained in Eqs. (E} and dHJ. 

At any event p on the spacetime manifold, we have a tangent space T p of vectors at p 
as well as the dual space T* of one-forms at p with corresponding expressions of the metric 
tensor given respectively by Eqs. fl2]) and fl9]). At p, the symmetric matrix {g^ u ) can be 
diagonalized, so that {g^ u ) = Q~ 1 DQ, where Q, QQ T = 1, is an orthogonal matrix and 
D = di&g(d ,di,d,2,d 3 ), where d$ < and di > 0, for i = 1,2,3. Thus by means of a 
local transformation and a certain scaling, (<7 M „) at each point p can be reduced to (rj^). 
Furthermore, we have that at p, (g^) = Q~ 1 D~ 1 Q, where the diagonal elements of D^ 1 
respect the metric signature, as before. In a similar way, (g^) can be reduced to (r] fiu ) at p. 
This pointwise duality thus stays at the threading level and does not extend to the slicing 
approach. That is, the conditions given in Eq. ([8]) are in general different from those given 
in Eq. ([6]), as can be seen via explicit examples discussed in Section IV. 

III. ADMISSIBILITY CONDITIONS 

The coordinate conditions fl6]) that result from the threading approach can be naturally 
incorporated into the underlying geometric framework of general relativity; however, the 
slicing conditions (|E|) can be neither invariantly formulated nor are they sufficiently local to 
become part of the structure of Lorentzian geometry. Indeed, all coordinate patches in which 
the metric takes the form g^ u {x) = r/^ + h^ v (x), where h^ u (x) can be treated as a small 
perturbation — such as in the Riemann normal coordinates about any event or the Fermi 
normal coordinates about the world line of an observer — satisfy both conditions ()6]) and OH]) . 
While coordinate systems that are compatible with conditions (|E|) protect chronology, this 
does not mean that chronology violation does not occur in the gravitational field under 



9 



consideration. In fact, general relativity per se does not require the existence of a cosmic 
time, as conditions (jSJ) cannot be woven into the fabric of Lorentzian geometry J4|. From 
this viewpoint, general relativity must be abandoned if one insists on integrating the slicing 
coordinate conditions into the foundations of the theory of gravitation. 

The threading and slicing approaches and their respective restrictions on coordinate 
charts have been discussed in one form or another by a number of authors. We should men- 
tion in particular the work of Zelmanov p, |6|, Cattaneo [7|, M0ller [8| and Lichnerowicz [9|, 
among others. A historical review is beyond the scope of this paper. Following Synge [101 ]. 
we refer to the combined threading and slicing conditions as the Lichnerowicz admissibilit 

n r 

requirements (see the theorem on page 9 of Ref. [9]). The book of Landau and Lifshitz [2 
contains a particularly clear discussion of the standard admissibility requirements. 

It is a theorem of linear algebra 11( that a real symmetric n x n matrix A is positive 
definite if and only if the principal minors of A are all positive. These are n scalars defined 

by 

"" A n ••• A lk 



det 



1 n. 



(17) 



Aki ■ ■ ■ Akk 

Now consider the standard threading admissibility conditions. These involve the matrix 
(9fj.u)' 9u < and 7^ = g^ — gugtj/gtt must form a positive definite matrix. Therefore, 



(7, 



711 


712 


713 


712 


722 


723 


713 


732 


733 



(18) 



must be such that 



7n > , 7n722 - 7i 2 2 > > det (7^) > . 



(19) 



Note that g u < and 711 = gn — g\ x j ' g tt > together imply that 



9tt 7n = det 



9tt gn. 
9it 9n 



<0 



(20) 
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Next, 7n722 — 7i2 > 0, when written out in detail, can be expressed as 



det 



7n 7i2 

721 722 



det 



9tt 



9tt 


9t\ 


9t2 


git 


9n 


912 


92t 


921 


922 



(21) 



Similarly, one finds that (see Appendix A) 

det(7»j) = 



9tt 



det(# 



hi/ j 



(22) 



Putting all these results together, we find that the standard threading admissibility condi- 
tions are equivalent to the statement that the principal minors of (g^ u ) must all be negative 
for our (—,+,+,+) signature. That is (gnu) must be a negative definite matrix. Simi- 
larly, the slicing admissibility conditions imply that (g^ u ) must be a negative definite matrix. 
The Lichnerowicz admissibility conditions then require that both (g^) and (g^ v ) be negative 
definite matrices, i.e., all their principal minors must be negative. 

In many spacetimes of interest in general relativity, Lichnerowicz conditions reduce to 
the standard admissibility conditions; indeed, most solutions of Einstein's equations are of 
this type 12] . This is the case, for instance, for the Kerr metric in the standard Boyer- 
Lindquist coordinates outside the static limit] however, CTCs occur in the region interior 
to the inner event horizon. Moreover, geodesic coordinate systems established along the 
world line of an accelerated observer in Minkowski spacetime are such that the threading 
and slicing conditions coincide 13|, Il4j; somewhat similar results are expected for radar 
coordinates |3|. The important point here is that Lichnerowicz conditions exclude CTCs, 
whereas the standard coordinate conditions do not in general exclude them. Much has been 
written about CTCs and time travel in general relativity — see, for example, 
recent review 19[ and the references cited therein. 



15 



18|, the 



IV. CAUSALITY AND WAVE PROPAGATION 

Consider the propagation of test radiation in a coordinate patch containing a CTC. The 
CTC might leave an imprint on waves propagating in such a spacetime domain and such 
a trace could then be observationally detectable. In this connection, it appears useful to 
study wave phenomena in spacetimes with CTCs. Consider, for example, a massless field 



11 

propagating in a spacetime region with a closed null geodesic. In the WKB approximation, 
we can imagine a wave packet moving essentially along the closed null geodesic, since in the 
WKB limit, massless wave propagation reduces to null geodesic motion. An examination 
of this situation reveals that the wave would appear to get immediately reflected: if the 
wave is initially moving forward in time, its spatial direction of propagation gets reversed 
when it starts to move backward in time. A similar result is expected for massive fields. We 
recall that a wave with propagation vector k^ = (u, k) in an inertial frame of reference is, 
up to a certain amplitude function, of the form exp(— iut + ik • x), so that as t increases, x 
increases along the wave vector k as the wave travels forward; however, if t decreases, the 
wave travels in the — k direction. When the direction of the arrow of time reverses, the wave 
appears to suffer a reflection. Thus we expect some unusual features for wave phenomena 
in gravitational fields that do not have a global time coordinate. 



Consider, for instance, the propagation of waves in the exterior Kerr spacetime [201 ]; 
in fact, there is absorption of radiation through the horizon in this case, but the scattering 
process contains little information about the CTCs that exist in the interior Kerr spacetime. 
Thus we must concentrate on wave propagation in regions where the CTCs are accessible. 



An ideal case would be the Godel universe 



2lM24 



24j , which is free of horizons and singularities. 



Do waves propagating in the Godel universe exhibit such an imprint [25H28J] ? In fact, the 
elucidation of this issue originally motivated the preliminary analysis presented in this paper 
and, as will become clear below, its complete treatment remains a task for the future. 

The metric of Godel's universe in the original quasi-Cartesian coordinates can be written 
as 

_ ds 2 = -dt 2 - 2V2Udtdy + dx 2 - U 2 dy 2 + dz 2 , U = e^* , (23) 

where Q > is the frequency of universal rotation and d t , d y and d z are three of the five 
Killing vectors of this spacetime. The determinant of the metric turns out to be g = —U 2 
and the inverse metric is thus given by 

- (d s ) 2 = (d t ) 2 - 2V2U-\d t )(d y ) + (d x ) 2 + U- 2 {d y f + (d z ) 2 . (24) 

The fundamental observers follow geodesies of metric ( l23l) : in fact, their four- velocity u is 
given by u = d t , so that we have 

G aP = tt 2 (g aj3 + 2u a u p ) . (25) 



12 

The source of Godel's solution of the Einstein field equations is dust of constant energy den- 
sity p, with fluid unit four-velocity aligned with the time coordinate lines, and cosmological 
constant A = —Q 2 = —Airp. Here Q > describes an intrinsic counterclockwise rotation of 
the universe around the z axis. The threading conditions are satisfied by metric ( 123|) . while 
it is clear from the inverse metric ( 1241) that g tt = 1. Indeed, the slicing coordinate conditions 
do not hold and the original Godel coordinates are therefore not admissible in the sense of 
Lichnerowicz. 

In 29[, exact Fermi normal coordinates were first explicitly introduced along the world 
line of a fiducial fundamental observer in Godel spacetime. Imposing Lichnerowicz conditions 
on these Fermi coordinates limits their domain of applicability to the interior of a cylindrical 



region of radius v21n(l + y/2)/Q ~ 1.25/f2 about the axis of rotation 29]. This boundary 
cylinder consists of closed null curves beyond which CTCs exist. 

Propagation of test electromagnetic radiation in the Godel universe was first investigated 



in [25|] with the express purpose of demonstrating the coupling of photon spin with the 
essentially uniform gravitomagnetic field of this rotating universe. The formal simplicity of 
the Godel metric made it possible to find the exact solution of the perturbation equations 
and thus establish the gravitational coupling of the helicity of the photon to the rotation 
of the Godel universe. To solve Maxwell's equations, it was deemed convenient to use the 
Skrotskii method, which is based on the analogy between a gravitational field and an op- 
tical medium 30], and replace the Godel field with an effective gyrotopic medium. As is 
well known, in the treatment of electrodynamics in arbitrary spacetime coordinates, it is al- 
ways possible to introduce instead a certain hypothetical medium with definite constitutive 



properties that occupies Minkowski spacetime in quasi- Cartesian coordinates 3JJ; in fact, 



this method is conceptually related to Weyl's formu 



relativity and to Gordon's optical metric [lOj, |32|, |33( . The physical description of photon 



ation of Fermat's principle in general 



helicity states is straightforward in scattering situations involving asymptotically flat space- 
times, but requires special care in universe models that are not asymptotically flat. For 
instance, in the Godel universe this connection is established via the limit as Q — > and 
Godel spacetime reduces to Minkowski spacetime. Furthermore, we expect that the spec- 
trum of electromagnetic wave perturbations of the Godel universe would contain certain 
discreteness properties in close analogy with the Landau spectrum of charged particles in 
a uniform magnetic field. Out of the complete spectrum of finite perturbations, there is a 



13 
part that satisfies these various physical requirements and only this part of the full spectrum 



was given in 25[ without providing details of the selection process. The preferred part of 
the spectrum involves wave vectors whose magnitudes are sufficiently large compared to Q, 
corresponding to a pulse of radiation produced by a localized source. From this solution one 
finds, in addition to the spin-rotation-gravity coupling, that waves cannot propagate parallel 
to the direction of rotation of this universe and, moreover, that waves can propagate only 
in one direction along the y axis. We wish to determine how these features as well as the 
spectrum are connected with the violation of causality in the Godel universe. 

Electroma gnet ic fields in the Godel universe have been discussed via the Debye potential 
formalism in 34j. Moreover, the scalar and neutrino perturbations of the Godel universe 



have been studied by a number of authors — see [35H39J] and the references therein. 

To see the imprint of acausality on these wave results, it is necessary to compare them 
with the results of wave propagation in causal Godel-type rotating universes that possess 
cosmic time. Indeed, Godel-type solutions of Einstein's equations exist that contain both 



causal and acausal universe models — see 40M43| and the references cited therein. The metric 



of such Godel-type solutions is of the form 

- ds 2 = -dt 2 - 2nR{t)Udtdy + R 2 (t)[dx 2 - {ri 2 - l)U 2 dy 2 + dz 2 ] , (26) 

where R(t) > is the scale factor, 

U = e^, V = 2 f. (27) 

Here A > is a constant parameter and Q > is a vorticity parameter. The Godel solution 
is recovered for R — 1, A = \/2Vt and i] = \/2. In general, Eq. (1261) represents a spacetime 
of Petrov type D with three Killing vector fields given by d x — Xyd y , d y and d z , so that it is 
a spatially homogeneous universe of type III in the Bianchi classification. Let 6^ = S^q be 
the four-velocity vector of the family of fundamental observers in the Godel-type universe, 
then this congruence has expansion and rotation, but no shear. In particular, the rotation 
tensor can be expressed as 

fa/3 = Q[a;0\ + Q>[ofip] , (28) 

where a M is the acceleration of this congruence given by a b = —r/RUdy with R = dR/dt. It 
follows that the only nonzero components of f a p are given by f xy = —f yx = £IRU, so that 



14 



the vorticity of the congruence of the fundamental observers can be expressed as 

n 



n 



nJapJ 



R{t) 



(29) 



Only positive square roots are considered throughout this paper. Moreover, the expansion 
is given by 3R/R. In the special case of stationary Godel-type universes with R — 1, the 
fundamental observers follow geodesies and the vorticity of their congruence is simply Q. 



Various sources have been considered for universes of Godel type 40 



431 ]; moreover, the 



corresponding two- fluid cosmological models have been discussed in [441 ] . 

For metric (1261) . we find g = —R 6 U 2 and hence the inverse metric is given by 

- (d s ) 2 = ( V 2 - l)(d t ) 2 - 2r ] R- 1 U- 1 (d t )(d y ) + R- 2 [(d x ) 2 + U- 2 {d y f + (d z ) 2 



(30) 



The threading coordinate conditions are satisfied here, while the slicing coordinate conditions 
are only satisfied for r\ < 1. Thus the solutions are causal for r\ < 1 and acausal for n > 1, 
while n = 1 constitutes a rather special limiting case. 

A detailed treatment of the propagation of electromagnetic waves in this Godel-type 
universe is due to Korotky and Obukhov [26] and Saibatalov |27j . who provided a complete 
description of the spectrum and outlined the main differences between the spectra in the 
causal a nd acausal cases. We emphasize that the more extended treatments of the Godel 
case in 



261 ] and |27j] go beyond the results of 



251 ]. where only the part of the spectrum 



that was deemed most phy sically reasonable was presented. In particular, the discrepancy 
between the results of [25] and |2j| is not due to any inherent limitation of the Skrotskii 
method. 

It is interesting to note that scalar and neutri no pert urbations of Godel-type universes 
have been the subject of recent investigations 



28 



45 



4&j. 



In the following section, we investigate the propagation of radiation parallel to the di- 
rection of rotation of a stationary Godel-type universe with R = 1. This two-parameter 
subclass has been thoroughly investigated in 47|, |48| and the references cited therein. 



V. WAVES IN STATIONARY GODEL-TYPE UNIVERSES 



Consider the propagation of waves in metric (1261) with R — 1. It turns out that the main 
result of this analysis is essentially independent of any particular wave equation; therefore, 
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we consider the scalar wave equation for the sake of simplicity. We seek a solution of 

V a d a ^ - mgtf = (31) 

in the form 

y = e i(-ut+k 2 y+k s z)^ x ^ ^ ( 32 ) 

where rfiQ = mo/h, mo is the mass of the scalar particle and u, &2 and k^ are constant 
propagation parameters. The function ip{x) satisfies 

U ~T i U lh ~ {*% + 27]Ujk2 U + ^ ~ 1)a;2 + k * + ™ol u2 }^ = ° • ( 33 ) 

Now introducing U^ 1 = x as a new variable and imposing the requirement that the wave 
travel in the z direction (fe = 0), we find 



v2 , 2 dy 

dx 2 



*V -T-, - [(^ " IV 2 + kl + m^ = . (34) 



Alternatively, Eq. ( 1331) with &2 = can be written in terms of variable x as a second order 
ordinary differential equation with constant coefficients and its solutions are of the form 
exp(— Xcrx). For x : — oo — > oo, x '■ °° — *" an d the solutions for ^ can be written as x° ■, 
where 

AV(ct-I) = {rj 2 - l)uo 2 + kl + ml . (35) 

The only finite solution for ip in the range x : — > oo is a constant, provided 

{rj 2 - l)u 2 + kl + m 2 = 0. (36) 

Therefore, if r\ < 1, so that the stationary Godel-type universe is causal, we have 

tf = ^ e *H"W-*a*) w 2 = M±^l ^ (37) 

1 — rj 

where ^o is a constant; otherwise, for 77 > 1 wave propagation parallel to the rotation axis 
is impossible. That is, waves can only propagate parallel to the rotation axis of causal 
stationary Godel-type universes. Let us recall here that waves can freely propagate parallel 



to the rotation axis of a compact gravitating source as well as a Kerr black hole [20|, |25|, 130 ] . 

What happens to a packet of scalar radiation produced in the neighborhood of z = zq in 

the acausal (rj > 1) Godel-type universe? It follows from Eq. (|36|) that the radiation is then 
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confined in space around Zq. That is, with k 3 = ±m, where we use the upper sign for z > z 
and the lower sign for z < zq, k > is given by 

(rj 2 - l)u 2 + ml] . (38) 

Then the solution for the scalar field can be expressed as 

q, = fa e - iut ~ K (*- zo) , z > z , (39) 



and 



* = Qe -i"t+K(z-z o ) ^ z <z , (40) 



where 0o is a constant. 

Consider now the propagation of radiation in the limit of high frequencies and wave 
numbers. For electromagnetic radiation, the geometric optics (or eikonal) approximation 
scheme results in the massless Hamilton- Jacobi equation and the corresponding propagation 
along null geodesic rays (49]. This is indeed the case for any reasonable massless wave 
equation. In the massive case, the analogous treatment is the WKB approximation, where 
for the scalar field we set 

f = e iS/R (*o + Wi + ...). (41) 

Substitution of this ansatz into Eq. ( l3i~|) results in an expansion in increasing powers of h 
starting with h~ 2 ; for H — > 0, the terms in the expansion proportional to h~ 2 , h~ l and h° 
must vanish. The primary result, which follows from setting the coefficient of h~ 2 equal to 
zero, is indeed independent of the nature of the wave and is the Hamilton- Jacobi equation 

g" v V IM SV v S + rr^ = 0, (42) 

where the scalar field S is related to action. The other equations then describe, for in- 
stance, how the wave amplitude \l/o propagates along the geodesic that is the solution of the 



Hamilton- Jacobi equation, and so on 49]. 



In the case under consideration, Eq. ( 1301) implies that 

For propagation along the z direction, it is natural to assume that 

S = -Et + P 3 z , (44) 
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where E and P3 are real constants representing the energy and momentum of the particle. 
Substitution of Eq. flU} into Eq. (g3]) yields 



( v 2 -l)E 2 + Pi + ml = 0. (45) 

This equation has no real solution in the acausal case (77 > 1), but in the causal case 
(77 < 1), the solution corresponds to Eq. ( 1371) . as expected. That is, the wavefront is the 
hypersurface of constant phase, which is proportional to S in the eikonal approximation; 
therefore, hypersurfaces of constant 5" must be spacelike (for massive fields) or null (for 
massless fields) in accordance with the Hamilton- Jacobi equation. 

In general, the solution of the Hamilton- Jacobi equation corresponds to a geodesic. To see 
this, let m w^ = V^S and note that w^ v = u>„ ;M ; that is, w^ is a curl-free vector field. Taking 
covariant derivative of Eq. ( )42l) . we find that g^ v w^ a w v = 0. This can be written, using the 
curl-free property, as g^ u w a ^w v = 0, which is the geodesic equation for w^ = dx^/ds. Thus 
in the case under consideration, m w^ = (-E, 0, 0, P3) and one can easily find the curl-free 
geodesic congruence in the causal case (77 < 1), namely, 

t = t + ^(l-r/ 2 )(r-r ), (46) 

where r is the proper time along the world line and 

x = x , y = y J i 77r^( r_r °)> z = z + —(t-t ). (47) 

m U(x ) m 

Here to, xo, yo, zq and To are constants of integration and E and P3 are connected via Eq. fj45l) . 
There is no corresponding geodesic for 77 > 1, as the Hamilton- Jacobi equation results in an 
imaginary energy in this case. 

It is puzzling that the geodesic corresponding to the eikonal limit of waves propagating 
in the z direction for 77 < 1 involves motion in the y direction as well. Let us note that in 
general, metric (I26"j) allows geodesies moving parallel to the axis of rotation regardless of the 
value ofrj; indeed, 

C = 7 (l,0,0,/3), (48) 

where (3 and 7 = (1 — /3 2 ) -1 / 2 are constants, is the unit four-velocity of such geodesic 
observers. However, m C/i is not the gradient of a scalar, since £ M is not a curl-free vector 
field, as can be directly verified. A bundle of neighboring geodesies has in general nonzero 
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vorticity; on the other hand, if the vorticity vanishes at one point, then the bundle will be 



completely free of vorticity — see page 84 of [23 ] . 



The remaining difficulty has to do with the fact that the Hamilton- Jacobi equation does 
not result in a real particle energy for rj > 1. This circumstance is a consequence of the 
violation of slicing coordinate conditions in metric ( 126]) . To see this in general, note that 
the geodesic equation is derived from the extremum of the action, 5S = 0, where 

S = -rn / ds (49) 

is a functional of the path that connects the two fixed events that are endpoints of the integral 
in Eq. (j4*9~l) . Let us define the Lagrangian L(t, x\ v^) = —mods/dt, where v l = dx % jdt. Using 
Pi = dL/dv l , the corresponding Hamiltonian H(t,x\pj) = PiV 1 — L works out to be — p , 
where po, as a function of t, x % and pj, must be obtained from g^p^Pv + ml = 0. This 
equation has the solution 

- Po ^( ( ,^ R 0^ P , ± (^±|^) 1/2 . (50) 



Thus the slicing coordinate conditions would ensure that the Hamiltonian in Eq. (150]) is 
real [50]. Then, integrating along geodesies in Eq. ( 149]) . but now with only the initial event 
kept fixed and the final event considered variable, we have 

S(t,x l ) = / (pjv 3 -H)dt, (51) 

so that dS = pidx 1 + podt, or V^S = p^, which then leads directly to the Hamilton- Jacobi 
equation. However, if Eq. (150]) cannot produce a real solution, then a real Hamiltonian does 
not exist and the connection of the geodesic equation with the Hamilton- Jacobi equation 
and hence wave propagation is severed. We should mention in passing that this derivation 
of the Hamilton- Jacobi equation implies that a fountain of neighboring geodesies emanating 
from a single event has vanishing vorticity. 

It is interesting to illustrate these results for the propagation of Dirac particles parallel 
to the rotation axis of stationary Godel-type universes. 

VI. DIRAC EQUATION IN STATIONARY GODEL-TYPE UNIVERSES 

Let us now consider the massive Dirac equation in a stationary Godel-type universe 
described by metric ( 126]) with R — 1. Here we follow the standard notation for the Newman- 
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Penrose formalism as employed in the monograph of Chandrasekhar 51] . Therefore, due 
to the nature of the subject matter, the notations and conventions used in this section 
are generally independent of the rest of the paper, except when otherwise indicated; in 
particular, we switch the sign of the metric signature. 

An orthonormal frame naturally adapted to the spacetime coordinates is given by 



ij = dt + r\Udy , u x = dx , ufi = Udy , u z 



dz, 



which can be used to form a null tetrad frame 



/ 



(or + u z 



n 



[us — U! 



V2" " V2 

The associated non-vanishing spin coefficients are 

V2 



m 



V2 



u + iuj 



u\ 



a = —(3 = 
while the only nonvanishing Weyl scalar is 



V2 



A, a — p — 2e = 27 = — 
4 ' P H ' 2 



n. 



V ; 2 



' 6 



(v 2 



(52) 



(53) 



(54) 



(55) 



as metric ( |26l) is of Petrov type D. Following 51], the massive Dirac equation in the Newman 



Penrose formalism is summarized by the following set of equations 

(D + e- p)F 1 + (5* + 7T - a)F 2 = im* G x , 

(A + fi- 7)^2 + (6 + P- r)F 1 = im* G 2 , 

(D + e* - p*)G 2 -(S + tt*- a*)d = im* F 2 , 

(A + p*- r f *)G 1 - (5* + /3*- r*)G 2 = im* F 1 , 



(56) 



where F 1; F 2 , G\ and G 2 are the components of the Dirac spinor, m* = m /\/2 is proportional 
to the mass of the Dirac particle (see [5l|, p. 543) and 



D = TV M , A = ^V M , 



S = m M V, 



(57) 



are directional derivatives along the null tetrad frame. Thus, 

(D + e - p)F 1 + (S* - a)F 2 = im* G x , 

(A + p - e)F 2 + (S- a)F x = im* G 2 , 

(D-e + p)G 2 -(5- a)d = im* F 2 , 

(A-p + e)G 1 - (6* - a)G 2 = im* F t , 



(58) 
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so that the spin coefficients are involved only through a real a and the combination 

iy/2 t 



e-p=-e 



-n. 



(59) 



which is purely imaginary Hence, Dirac equation can be finally expressed as 



(D - e)F x + {6* - a)F 2 

(A + e)F 2 + (6- a)F 1 

{D + e)G 2 -{8-a)G 1 

(A - e)Gx - {5* - a)G 2 

Let us look for separable solutions of the form 



im* G x , 
im* G 2) 
im* F 2 , 

im* F\ . 



(60) 



[Fi,F 2 ,Gi,G 2 ] = [X l (x),X 2 (x),Y 1 (x),Y 2 {x)]e 



-iujt J rik2y+ikzz 



(61) 



with k 2 = 0, as for scalar waves examined in the previous section. 

The resulting equations for the spinor components X\ y2 and Y\ >2 form a first order system 
of ordinary differential equations with constant (complex) coefficients, namely, 

( X x \ I A + B + b \ { Xx\ 
B_ A_ b 

-b A + C + 
\-b C_ A_ ) 



d 
dx 



X 2 
Yi 



X 2 
Y-i 



(62) 



where b = irh and 

A± = --±ur), B ± 



-u 



± h + - 



C 






to 



i ft 
± fo - - 



(63) 



Let us note that the function U = exp(Ax) appears in Eq. fl52l) and the subsequent analysis 
in association with variation along the y direction, but such variation disappears once we 
set k 2 = in Eq. (IBID. 

The eigenvalues of the 4x4 matrix in Eq. ( J62l) determine the solutions. These are given 

by 

~u\r? - 1) H 



Ai. 



2 



mg + kl + E'^ 



,/ft 
2 



1/2 



(64) 



where E and E' independently can be either +1 or — 1; that is, E 2 = E /2 = 1. The elementary 
solutions of system (162]) are of the form e AlX , where I — 1, . . . , 4. We require that the Dirac 
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equation have finite solutions for all the spinor components for x : — oo — > oo. The only 
possibility is then that of having constant solutions for X 12 and Y lf2 , i.e., the determinant 
of the 4x4 matrix in Eq. ( 1621) must vanish. Thus the allowed values of u satisfy the relation 

2 \2 



"V-i)+(v^+*i±i) -t = ' (65) 



which agrees with Eq. (136]) in the WKB regime. That is, with E = hu and P3 = M3, 
Eq. (|65|) can be written as 

(r] 2 - 1)E 2 + Pi + m 2 ± m(Pl + ml) 1 ' 2 + h 2 n 2 (l - A^ 2 ) = , (66) 

so that as h~ — > 0, while E, P3 and rriQ are held fixed, we recover the WKB limit. 

The two branches of Eq. (|66|) reflect the spin 1/2 character of the Dirac particle; indeed, 
with m = 0, the helicity-rotation coupling evident in the spectrum ( 166]) is precisely reminis- 



cent of the corresponding result for photons 25| . when one allows for the different helicities 
involved. 

For a localized packet of cosmic ray protons, for instance, we expect that the packet is 
dominated by k 3 » A and k^ » Q, where 1/A and 1/Q are length scales associated with 
the Godel-type universe. Thus propagation of the packet parallel to the rotation axis is 
generally possible in the causal case (77 < 1) and impossible in the acausal case (77 > 1). In 
the latter situation, we expect confinement of the high- energy particles in space, just as for 
the scalar field — see Eqs. (15§|) and (14"U|) — since waves with /c 3 < A and k 3 < Q may be able 
to propagate in accordance with Eq. ( IS"5lh 

VII. DISCUSSION 

The most basic measurements of an observer consist of the determination of temporal 
and spatial intervals. In this connection, we have explored the admissibility of coordinate 
patches in some open spacetime region on the basis of splitting spacetime into its elements 
via the threading and slicing approaches. All coordinate systems in general relativity must 
satisfy the threading (Landau-Lifshitz) conditions. The Lichnerowicz admissibility condi- 
tions involve the additional slicing coordinate requirements, which protect chronology and 
exclude CTCs. 

If it were possible to incorporate the slicing coordinate conditions into the geometric struc- 
ture of Lorentzian geometry, violations of chronology and hence CTCs would be forbidden 
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in the general theory of relativity. However, since the "local" t = constant hypersurfaces 
cannot be naturally absorbed into the foundations of general relativity, the slicing conditions 
simply place additional restrictions on possible coordinate systems. 

The threading and slicing coordinate requirements together ensure that the coordinates 
assigned to spacetime events are physically reasonable. As such, coordinate systems that 
can be physically constructed, such as the GPS system, are expected to satisfy Lichnerowicz 
coordinate conditions. We have examined the physical implications of Lichnerowicz condi- 
tions by studying wave propagation in spacetimes that do not possess universal temporal 
coordinates. In such stationary rotating Godel-type universes, for instance, we show that 
radiation in the WKB approximation cannot propagate parallel to the axis of rotation. Pos- 
sible generalizations of this result as well as further consequences of chronology violation for 
wave propagation are topics for future investigation. 
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Appendix A: General Form of Equation (j22j) 



Imagine an arbitrary (1 + n) x (1 + n) matrix of the form 



^4oo An • • • Ao n 



A 



A^IO ^11 • • • A\ n 



A n o A n \ ... A r , 



(Al) 



which can be denoted as 



.4 



Aio Aij 



i,j = 1,2,..., n. 



(A2) 



We assume that Aoo 7^ 0; in fact, this condition can be ensured — unless A is trivial in the 
sense that it has only zeros in its first row and column — by appropriately shifting rows or 
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columns, as we are only interested in det(A). Then, it is possible to show that 

det(A) = A 00 det (Aj - ^i) . (A3) 

Consider the first row of matrix A, 

R = [A 00 , A 01 , . . . , A 0n ] . (A4) 

Now subtract R (A 10 /A Q0 ) from the second row of matrix A, R (A 2 o/A 00 ) from the third 
row and so on. In this process, the determinant of the matrix does not change due to its 
alternating character but matrix A is transformed into 

B=\ A ' a *"._... I. (A5) 



A 



A m A, 



"j 



% 1 A 00 

The determinant of this matrix is simply 



dci-(B) = A l){) dvt{A ij -^ i ) ■ (AG) 



which proves Eq. (1A3|) . since det (A) = det(S). 
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